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ABSTRACT 


Perturbation theory for fixed sources was applied to radiation shielding pro- 
blems to determine changes in neutron and gamma ray doses due to changes in 
various shield layers. For a given source and detector position the perturba- 
tion method enables dose derivatives due to all layer changes to be determined 
from one forward and one inhomogeneous adjoint calculation. The direct approach 
requires two forward calculations for the derivative due to a single layer 
change. Hence, the perturbation method for obtaining dose derivatives permits 
an appreciable savings in computation for a multilayered shield. For an illus- 
trative problem, a comparison was made of the fractional change in the dose 
per unit change in the thickness of each shield layer as calculated by pertur- 
bation theory and by successive direct calculations; excellent agreement was 
obtained between the two methods. f 


INTRODUCTION 


Reactors for space propulsion systems or auxi 1 i ary power systems require 
shielding to protect personnel, system components, or payload from neutrons 
and gamma radiation. These shields, are usual ly„composed_ of alternate layers 
of heavy metal and hydrogeneous material. The shielding problem for these 
systems is to determine the number, arrangement, and thicknesses of these shield 
layers that will result in minimum weight shields while maintaining the radiation 
field at various positions around the system within prescribed limits. 

One-dimensional shield weight optimization procedures have been evolved for the 
design of layered configurations (refs. 1, 2, 3) and these have been extended 
to constraints in several directions around two-dimensional configurations 
(refs. 3 } 4, 5). All of these require the determination of the change in dose 
at a detector due to a change in thickness of each layer of a given shield con- 
figuration. These are generally obtained with Sjj transport calculations requiring 
at least as many separate calculations as there are shield layers. As the op- 
timization proceeds and the configuration changes, this procedure has to be re- 
peated unt i 1 the conf i guration converges to an optimum one. The optimization 
procedure becomes time consuming and particularly for the two-dimensional case 
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would require excessive amounts of computer time. 

The present work applies perturbation theory for systems with fixed sources 
to this shielding problem. The perturbation method utilizes only one forward 
and one inhomogeneous adjoint solution for an initially selected layered shield 
configuration to determine all necessary dose derivatives at a detector for a 
given neutron or gamma ray source component. The method appreciably reduces the 
amount of calculation required for multilayered shields. In addition, pertur- 
bation theory yields a considerable amount of information concerning the effect 
of various physical processes (elastic scatter, inelastic scatter, etc.,). on 
the dose rate. 

The perturbation method in this work is applied to a simple spherical system 
to determine the dose derivatives of the various layers. This illustrative 
problem also serves to demonstrate the ability of the perturbation method to 
extract additional informative data concerning interaction processes of the 
various shield materials. 


ANALYSIS 


The theory underlying the application of perturbation methods to non-reactivity 
problems appears mainly in the Russion nuclear energy literature (refs. 6 and 
7). Lewins reviews in reference 8 developments in perturbation theory. This 
section presents the equations used in the perturbation analysis for systems 
with fixed sources. The theory will be outlined for a fixed source neutron 
system since the results for fixed source gamma ray system are similar in form. 


Consider the time-independent multigroup form of the inhomogeneous neutron trans 
port equation. This equation is given by 





* v ♦ r » 


a. 


(i) 


The notation used for equation (1) is conventional: ) is the angular 

flux £or group g, ^ is the macroscopic total cross section for group g, 
Qg(r^H) is the angular fixed source distribution for group g,G^(g^g) is 
the X. -th order expansion coefficient of the scattering from group g' to group 
g, <&&*», GT) i s the ,£,m-th expansion coefficient of the;.angular flux in a spher- 
ical narmonics series, and is the ^m-th spherical harmonic as a 

function of the unit vector j£ . The boundary condition for the angular flux 
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at the convex boundary of the system is taken as 

" D i ( 2 ) 

This boundary condition is the usuaj condition implying that no neutrons enter 
the system from the vacuum. Here indicates the outer boundary of the system 
and JC des i gnates an inwardly directed direction 

Once the solution to equation (l) has been obtained, then a physical process 
which depends linearly on the angular flux may be written as a functional F. 

The multigroup form of F is given by 

P- 2^4^ OtjfC) ^ oljldlL ^ 

The functional F gives contributions caused by the interaction process for 
the neutrons in the system. The quantity f has the role of an inter-? 
action cross section for group g at position 7? for neutrons moving in the 
d i rection JL 

Following Usachev (ref. 6) an equation can be derived fqjr the importance of 
neutrons for the physical process described by the functional F. This equa- 
tion for the importance which is equal to the effect F at some 

detector j>osi t ion or region due to a unit source of neutrons placed in group 
g at ), is given in multi group form by 



r. • # t G / (4) 

The boundary condition for the importance follows from the fact that, 

if a neutron escapes from a system, it cannot contribute to the physical pro- 
cess of interest. Therefore its importance must be zero and the boundary- 
condi t i on becomes 

= O ; ’" ,Gr (5) 

where XlT designates an outwardly directed direction ,(L . 


4 


The next step in the analysis is to derive a perturbation equation in terms of 
a perturbation in the density of the material constituents of the system. The 
perturbation in the density of a material, which can be related to a perturbation 
in the thickness of a given shielding layer, is introduced in the equations through 
the macroscopic cross sections and appears as a perturbation in the cross sections 
of the system. The solution of equation (l) with the perturbed cross sections 
yields a perturbed value for the angular f 1 ux ^ ) . The equation for the per- 
turbed system is 



along with the boundary condition 

J<3 Gia, sr) -o ~ - - • £ 

In the above equation represents the perturbed value of the angular 

flux for group g. Perturbations in the cross sections are indicated by S(T. 

Also the source term Qg(/2,.£L ) is assumed to undergo a perturbation in the amount 
of § 0.9 ) . 

The neutron transport equation for the perturbed system, equation (la), along 
with the equation for the importance can be manipulated to give the fallowing 
result:- . 



- s fa fa * 04*0 + %<%(*, .*>] 

^ a ^ 

+ r. fa fa ^ = ° « 

a 
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Equation (6) is an exact relationship and represents a condition which the perturbed 
angular flux and importance must satisfy. If all of the perturbations in equation 
(6) are zero, then the equation reduces to 

f - S ^ A XL 
% Jt 

— S ^ fejil) Gfy ft; a) Qg.^ 

^ Ji. 

has been replaced by its unperturbed value 
expressed by equations ( 3 & 3a) is well known. 

- P 

In the perturbed systems the functional r giving the contributions for the 
physical process of interest is equal to 




The result 



^ J/£ ^ ft, ft) P^ftjft) 


a 

Here the general case of a perturbed system is considered for which fey i/L ) A 
is the perturbed value of the quantity faftyit) . The quantity f is 
taken to be 0 



ft) 





_ p ^ 

Next consider the difference af of the functionals p and Using equations (3a) 
(7) and (8) in equation (6) results in 



Equation (9) is the desired perturbation expression for a system with fixed 
sources. This equation gives the change in the functional F for a given physical 
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process due to changes in various neutron cross sections, changes in the fixed 
source, and changes in the physical process itself. Equation (9) is an exact 
relationship since no approximations are inherent in its derivation. However, 
since only first order perturbation theory is being considered, the perturbed 
flux may be replaced by its unperturbed value 


J_ji shielding analysis, generally the biological dose rate at a detector position 
KA is the functional of greatest interest. The interaction cross sections 
ffoOtjjt) are the flux to dose conversion factors which are independent of 
Thus, we can express tg as 


0 ’ M 


W 

where dg is the flux to dose conversion factor for group g and 

is the Dirac delta function. Insertion of equation (10) into equation (3) 

9 i ves 




p - f h ^ ^ 


- ” Qi) 


Thus the functional F is 6 flt-ff) times the dose rate at position Jljt , 

ILLUSTRATIVE PROBLEM 


The spherical system chosen for this study shown in figure 1 has a central 
cavity or void region with a radius of 30 centimeters. This void region is 
followed by a 15 centimeter thick layer of natural tungsten and a 60 centimeter 
thick layer of 1 i th i urn hydr i de. The dimensions and composition of this system 
are 1 i s ted in Table I . 


The neutron dose was calculated using 26 energy groups at a detector located 
0.25 centimeters from the outer boundary of the system. This group structure 
is shown in Table II. The fixed neutron source of 1.7 x 10^7 neutron/seq was 
uniformly distributed in the central void region. The neutron fluxes and im- 
portances were determined using the Sn method with elastic scattering treated 
through the P3 order and using an S ] 6 Gauss-Legendre angular quadrature. The 
effect on the neutron dose as measured at the detector position was then deter- 
mined for tungsten and lithium hydride for all the important neutron reactions 
by using the perturbation formula given by equation (9), The results obtained 
with the perturbation formula are for a nominal perturbation in the density of 
each layer. 

The total effect of the tungsten or lithium hydride layer can be expressed as 
the fractional change in the neutron dose functional per unit change in layer 
thickness. These quantities are compared with similar results obtained from 
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successive direct flux calculations in which a nominal change in the density 
(or equivalently thickness) of a layer is made. 

Neutron interactions with tungsten nuclei generate a secondary source of gamma 
rays throughout the tungsten layer. The secondary gamma ray sources are due 
to the W(n,Tf) and W(n, n 1 ^ ) reactions. Results for these two secondary gamma 
ray sources will be presented separately. The analysis and results obtained 
for these secondary gamma ray sources follows that presented for .the fixed 
source neutron problem; however, perturbations in the secondary source due to 
perturbations in the W(n,Tf ) and W(n,n'2f) reactions caused by perturbations in 
the W density must also be included. All of the gamma ray calculation were 
obtained using the group structure given in Table III and for an Si 6 Gauss- 
Legendre angular quadrature with scattering treated through the P 3 order. 


RESULTS 


Effect of the Tungsten and Lithium Hydride Layers on the Neutron Dose 

The neutron dose rate was calculated from equation ( 3 ) to be 2.90 x 10 ^ rem/hr 
for the source of 1.7 x 10^7 neutrons/sec. The neutron dose calculated using 
equation (3a) was computed to be 3-11 x 10^ rem/hr. The ratio of the dose cal- 
culated from equations (3a) and ( 3 ) was equal to 1 . 07 - This ratio should be 
unity and departure from unity is a measure of the approximations inherent in 
the Sn method for a given spatial mesh and angular quadrature. Thus, if the 
values of the dose functional differ by a large enough amount to preclude 
obtaining the accuracy desired, then S n calculated angular fluxes and importance 
must be recalculated with a finer spatial mesh and angular quadrature. For the 
illustrative problem being considered, the error in the perturbation results is 
approximately half the percent difference rn:the values of the two dose functionals. 

Table I V presents the effect on the^neutron dose functional, kV of the tungsten 
layer. The data in Table IV gives the effect on 477&-P of each reaction which 
a neutron can undergo with nuclei of tungsten. Inelastic neutron scattering 
accounts for 60.9 percent of the contribution to the total kTS-F; elastic neutron 
scattering contributes 19.6 percent; (n, 2 n) reactions contribute 18.8 percent; 
and neutron absorptions contributes 0.7 percent. The data in Table IV indicates, 
as: an example, that a 1 percent increase in the density of the tungsten layer 
would decrease the dose of 2.90 x 10 ^ rem/hr at the detector position by 10. 3 
rem/hr. Also given in Table IV is the change in the neutron dose functional 
per gram of material. A value of -2.00 x 10"^ rem/hr/gm was obtained for W. 

Table IV also gives the change in the neutron dose functional per centimeter. 

For natural tungsten this coefficient is -68.9 rem/hr/cm. 

Table V gives a comparison of the fractional change in the neutron dose func- 
tional per unit change in the W layer thickness,Vf/p&t between the perturbation 
theory result and the result calculated by successive direct transport calcula- 
tions. The perturbation theory result is -.238 cm - l while the result from direct 
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flux calculations is -.234 cm - '. The two results are in excellent agreement 
wi th each other. 

Table VI presents the effect on the neutron dose functional, 47fSF ;> of the 
lithium-hydride layer. In addition, partial effects due to constituents H 
and Li are also presented. »Table VI indicates the dominant role of elastic 
neutron scattering for low mass number elements on the dose functional. The 
data in Table VI shows that a 1 percent increase in the lithium-hydride density 
would decrease the dose of 2.90 x 10^ rem/hr at the detector position by 24.0 
rem/hr; hydrogen and lithium contribute 12.2 and 11.8 rem/hr, respectively, of 
the 24.0 rem/hr. Table VI also indicates that 1 gram of LiH, uniformly dis- 
tributed throughout the LiH layer, reduces the neutron dose by 7-19 x 10 -if 
rem/hr. Also shown in Table VI is the dose coefficient per centimeter for 
LiH. For the LiH layer this coefficient is -40.0 rem/hr/cm. 

Table V gives a comparison of the fractional change in the neutron dose func- 
tional per unit change in the LiH layer thickness, SP/FSt;, between the pertur- 
bation theory result and the result calculated by successive direct transport 
calculations. The perturbation theory result is -.138 cm - ' while the result 
from direct flux calculations is -.140. The two results are in excellent agree- 
ment with each other. 


Effect of the Tungsten and Lithium Hydride Layers on the Photon Dose 


The |hoton dose rate was calculated from equation (3) to be 2.59 x 10^ rad/hr 
for a total secondary gamma ray source, due to neutron absorptions in W, of 
3.98 x 10'7 photons/sec. Using equation (3a) the photon dose was also calculated 
to be 2.59 x lcA rad/hr. The photon dose rate was calculated from equation (3) 

to be 7.10 x 1C)3 rad/hr for a total secondary gamma ray source, due to neutron 

inelastic scattering in W, of 8.17 x 10'7 photons/sec. Using equation (3a) the 
photon dose for thLs _ca S e was., calculated to, be J. 19- x 1 03 rad/hr. . 

Table VII presents the effect on the photon dose functionals, 4^^ , of the 

tungsten layer. The data in Table VII gives the effect on 4'77SP of each reaction 
which a photon can undergo with a nuclei of W; in addition, the effect on 4?7 r Sr 
of the perturbations in the secondary photon source are also given. From the 
data in the table it can be seen that a major contribution to 4 is due to 

these perturbations in the secondary photon source. 

Table VIII gives a comparison of the fractional change in the photon dose func- 
tional per unit change in:the W. layer thickness, , between the perturbation 

theory result and the result calculated by successive direct transport calcula- 
tions. The data in the table shows that there is excellent agreement in the 
values calculated by the two methods. 

Table IX summarizes the results obtained for the effect of the LiH layer on the 
photon dose due to the two secondary gamma ray sources. From Table VIII it can 
be seen that the coefficient S-p/f&f computed by perturbation theory is in excellent 
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agreement with the result obtained from successive direct flux calculations. 

CONCLUSIONS 


In summary, the results obtained for the coefficient Sf/fSt by using perturba- 
tion theory for fixed sources were in excellent agreement with those obtained 
using successive direct transport calculations. Perturbation theory also gives 
a considerable amount of information concerning each reaction which a neutron 
or photon can undergo with the nuclei of various materials and thus add to the 
qualitative understanding of changes in shield design. This same information 
is not easily or readily obtainable from successive direct transport calculations. 
The extension of perturbation theory to two-d imens ional geometry would pose no 
additional difficulties other than those of the transport solution itself. 
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FIGURE 1. Sketch of spherical system used in illustrative problem. (The detector 
position is 0.25 centimeters from the boundary.) 
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TABLE I. Size and Composition of Spherical System 
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TABLE II.' Neutron Energy Group Boundaries, Normalized Fission Spectrum, and Neutron 
Flux-to-Dose Conversion Factors. 


GROUP 

UPPER ENERGY BOUNDARY, 
eV 

NORMALIZED FISSION 
SPECTRUM 

NEUTRON FLUX-TO-DOSE 
CONVERSION FACTOR 
(rem/hr) /(neuts/cm 2 /sec) 

1 

1.492 x 10? 

.000155 

1.5 x 10-4 

2 

1.221 " 

.000886 

1.5 

3 

1.000 " 

.003451 

1.5 

4 

8.19 x 10 6 

.018725 

1 .48 

3 

6.07 

.029349 

1.4 

6 

4.97 

.047774 

1.35 " 

7 

4.07 

.106381 

1.3 " 

8 

3.01 

.088247 

1.28 

9 

2.47 " 

.046962 

1.28 

10 

2.23 

.095373 

1.28 

1 1 

1.83 

.136033 

1.35 " 

12 

1.35 

.080244 

1 .4 " 

13 

1.11 

.069595 

1.35 " 

14 

9.07 x 105 

.127183 

1.0 

15 

5-50 

.049062 

7.6 x 10" 5 

16 

4.08 " 

.085216 

5.0 

17 

1.11 " 

.015364 

2.0 " 

18 

1.50 x 10 4 

0 

0 

19 

3.35 x 103 

0 

0 

20 

5.83 x 10 2 

0 

0 

21 

1.01 

0 

0 

22 

2.90 x 10 1 

0 

0 

23 

1.07 

0 

0 

24 

3.06 x 10° 

0 

0 

25 

1.13 

0 

0 

26 

4.14 x 10" 1 

0 

0 


MOTE: Lower energy boundary of thermal group is 0.001 eV. 
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TABLE ill. Photon Energy Group Structure, Secondary Photon Spectra for Tungsten, 
and Photon Flux-to-Dose Conversion Factors 



Upper Energy 
Boundary, eV 

Number of 
Photons/abs 
i n W 

Number of 
Photons/inel . 
event in W 

Photon Flux-to-Dose 
Conversion Factors 
( rad/h r^pho tons /cm^- sec) 

1 

8.00 x 10 6 

.1067 

.00030 

6.510 x 10 -6 

2 

5.50 

.0610 

.00038 

5.712 

3 

5.00 

.0547 

.00042 

5.345 

4 

4.50 

.0471 

.00047 

4.895 " 

5 

4.00 

.2000 

.00053 

4.513 

6 

3.50 

.2308 

.00154 

4.102 " 

7 

3.00 11 

.2071 

.00714 

3.697 

8 

2.60 

.2417 

.01250 

3.322 ", 

9 " 

2.20 " 

.3250 

.03000 

2.913 

10 

1.80 

.2857 

.44444 

2.472 " 

1 1 

1.35 

.4444 

.88889 

1.883 " 

12 

.90 " 

.7538 

.07692 

1.105 

13 

.40 

.3636 

.15152 

5.886 x 10-7 

14 

.26 

0 

0 

3.240 " 

15 

.15 

0 

0 

0 


NOTE: 


Lower energy group boundary of group 15 


is 8.0 x 10 4 eV. 







TABLE IV. Neutron Dose Functionals and Derivatives for theTTungsten Layer 


Change in Dose Functional for Each 
Neutron Reaction, 4 frSP , rem/hr 

Elastic scattering 
Inelastic Scattering 
(n,2n) Reaction 
Absorption 

- 2.03 x 10 2 
- 6.30 x 10 2 
-1.94 x 10 2 
- .07 x 10 2 

Total Change in Dose Functionals 
4?T&fS rem/hr 

-10.3 x 10 2 

Change in Dose Functional per Gram 

-Z.00 x 10" 4 

rem/hr/gm 


Change in Dose Functional per cm. 

-68.9 

rem/hr /cm 
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TABLE V. Comparison for the Neutron Problem of Deri vat i ve Calculated by 

Perturbation Theory and Successive Direct Transport Calculations. 


Fractiona I 
Increase in 

Dens i ty. 


W 


LiH 




,02 


,01 


Equivalent Thick- 
ness, 


.30 

.60 


Fractional Change in Neutron 
Dose for Unit Change in Layer 
Thickness, -> 


i n 
CsrrC^ 


Direct Calc. I Perturbation Calc 


,234 


-.140 


.238 


-.138 


Pert. Calc . 
Di rect Ca l a 


1.02 

.99 
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TABLE VI. Neutron Dose Functionals and Derivatives for the Lithium-Hydride Layer 


MATERIAL 

H 

Li 

Change in. Dose Functional for Each 
Neutron Reaction, rem/hr 

Elastic Scattering 
Inelastic Scattering 
Absorption 

Total Change in Dose Functional, 

4 If'St ' , rem/hr 

-12.2 x 10 2 

-6.99 x 10 2 

-4.70 x 10 2 
- .09 x 10 2 

-12.2 x 10 2 

-11.8 x 10 2 

Change in Dose Functional per Gram, 
rem/hr/gm 

- 2 . 90 x 10" 3 

-4.04 x 10 _i+ 

Change in Dose Functional per cm, 
rem/hr/cm 

-20.4 

-19.6 


















Secondary Photon Source 


Change in Dose Functional for Each 
Photon Reaction, klf S-? , rad/hr 

Scattering 

Absorption 


Change in Dose Functional for Second- 
ary Source Perturbations, kV SF~ , 
rad/hr 


Total Change in Dose Functional, 
4 rem/hr 


Change in Dose Functional per gram, 
fad/hr/gm 


Change in Dose Functional per cm, 
rad/hr/cm 


W(/n^ ) Reaction | W(/»1 jVy) Reactions 

-7.20 x 10 3 
-2.35 x 103 

-8.80 x 10 3 
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TABLE VIII. Comparison for Photon Problem of Deri vat i veSf/p&tCalcul ated by Per- 
turbation Theory and Successive Di rect Transport Calculations. 


)■ 



REGION 

Fract i ona 1 
Increase in 
Dens i ty. 

Equivalent Thick- 
ness, 

Fractional Change in Photon 
Dose per Unit Change in Layer 
Thickness, & F/F&t ) CwC 1 

Pert. Calc. 
D i rect Ca lc 

'1 


Di rect Calc. 

Perturb. Calc. 


\ 

Secondary 

Photon 

Source; 

W ) 

Reaction 

W 

.02 

• 30 

-.107 

-.110 

1.02 


liH 

.01 

.60 

-.026 

-.026 

1.00 


Secondary 

Photon 

Source; 

) 

Reaction 

W 

.02 

• 30 

-.168 

-.172 

1 .02 


LiH 

.01 

.60 

-.039 

-.139 

1 .00 
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TABLE IX. Photon Dose Functionals and Derivatives for the Lithium-Hydride Layer 


Secondary 
Photon Source 


Change in Dose 
Functional for 
•each Photon 
Reaction, 

, rad/ hr 

Scattering 

Absorption 


WfoOf.- )Reacti on 


-l.OIxlO 4 -3.02x10 4 
- . 77*1 0 2 -5. 95x1 0 2 


Total Change in 

Dose Functional, -1.02xl0 4 -3.08xl0 4 

4'OfS-F - , rad/hr 


Change in Dose „ ' y 

Functional per -2.41x10 -1.05x10" 


gram, rad/ hr^n 


Change in Dose 
Functional per 
cm, rad/hr/cm 


2 -5. 1 3xl0 2 


W (rt\) / *i ) Reaction 


H Li 


■4. 02x 1 0 4 

■6.72xl0 2 


-4.11x163 -1.23X10 4 


-16.2 


-4.09X10 4 -4.11x103 - 1 . 2 3x 1 0 4 


-1.22x10-2 -9. 75x10“ 3 -4.22x1 0” 3 


■6.82xl0 2 -68.49 - 2 . 05 x 1 0 2 


-4.92xl0“3 


-2. 74x1 0 2 



























